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ABSTRACT 

The aim of this paper is to compare the two cases of an isolated neutron star, with a 
non-accreted crust, and that of an accreting neutron star, with an accreted crust, and 
try to estimate which one of the two would make a better source of gravitational waves. 
In order to do this we must evaluate the maximum "mountain" that the crust can 
sustai n in these two cases. We first do this using the formalism of lUshomirskv et al.l 
and find that the maximum quadrupole is very similar in the two cases, with the 
non-accreted crust sustaining a slightly larger mountain. We then develop a perturba- 
tion formalism for the problem, that allows us to drop the Cowling approximation and 
have more control over the boundaries. The use of this formalism confirms that there is 
not much difference between the two cases, but leads to results app roximately one or- 
der of magnitude larger than those we obtain with the formalism of lUshomirskv et al.l 

El 



1 INTRODUCTION 



The possibility that the crust of a rotating neutron star may sustain a "mountain", thus crea t ing a source of gravitational 
wa ves, has recently spa rked much interest. T h is is m ainly due to the suggestion bv iBildsterJ il998l) . echoing earlier ideas 
by IWagoneJ (|l984h and IPapaloizou fc PringQ l)l978l) . that the gravitational wave torque due to such a mechanism may be 



dictating the spin e quilib rium period for the LMXBs. A recent estimate of how large a qu adrupole one may expect was given by 
lUshomirskv et all ll200(j) . using the perturbation formalism of lMcDermott et alj j 19881) . From their work it appears that the 
quadrupole may be big enough to balance the accretion torque, if the breaking strain of the crust (a very poorly constrained 
parameter) is taken to be near the ma^dmum of the possible values. Even th ough other gra vitational wave mechanisms have 
been sugges ted for these systems (toroidal magnetic fields were suggested bv ICutleJ l)20n2l) and r-modes bv lAndersson et"ail 
l)l999ll200ri )). and even though one can build models of the interaction betwe en the accretion disc and the magnetic field that 
give the correct distribution of spin frequencies, as in lAndersson et all 1)20051) . the matter is still of great interest. 

There are in fact many questions that can still be asked, the main one being if there is a difference between the "mountain" 
that one can build on the non-accreted crust of a rotating neutron star and the one that can be built on the accreted crust of a 
neutron star in a binary system (such as the LMXBs). We shall see that these two different scenarios can lead to significantly 
different compositions for the crust. The case of isolated neutron stars is particularly interesting a s current gravitatio nal wave 
detectors are beginning to put astrophysically relevant limits on the ellipticities of known pulsars ([Abbot et all (|200i4l ). In this 
work we will discuss this issue, and present a new scheme for calculating the maximum quadrupole which allows us to have 
better control of the boundary conditions at the base of the crust than in the scheme of lUshomirskv et al.l |2Qp7). Another 
important point is that we do not adopt the Cowling approximation, which, as we will show, can lead to significant differences 
in the results. We intend to present a formali sm in terms of per t urbat ions of a spherically symmetric crust, without making 
the Cowling approximation that was made in lUshomirskv et alJ (|200(t . 

We will consider two cases for the crust core interface: in one case we will solve the hydrostatic equilibrium equations in 
the crust and in the core, and impose continuity of the tractions at the junction; in the other case we consider the core to 
remain unperturbed, and thus spherical. This is not physical but, as we shall see, allows us to estimate the errors associated 
with neglecting the perturbations of the core. The main point, however, is that if we leave the core unperturbed, this allows us 
to solve the TOV equations to obtain it's structure; we can thus obtain an estimate of how important relativistic effects may 
be. Such effects are relevant because, even though the Newtonian approximation is excellent at the low densities of the crust, 
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the structure of the core can be highly relativistic and have an impact on global parameters such as the mass and radius of 
the star. We shall first develop the formalism for an n = 1 polytrope and then apply it to the case of a realistic equation of 
state, for an accreted and a non-accreted crust. 



2 MAXIMUM MOUNTAIN AND CONSTANT STRAIN IN THE CRUST 

In this section we shall briefly review the method of lUshomirskv et al.l J2000l) . who were the first to perform a detailed 
calculation to estimate the maximum size of a "mountain" that a neutron star crust could sustain. As it will be necessary to 
refer to it, let us recapitulate briefly their work. 

Let us assume that the crust responds elastically to pressure and density perturbations. We can thus treat it as a solid 
with a shear modulus /^(r) and write the stress tensor as: 

(1) 



Tab 



-pQab + tab 



where Qab is simply the flat 3-metric and tab is the shear stress tensor of the crust, which vanishes in the fluid interior. We 
will consider the equilibrium shape of the star to be spherical, thus effectively representing the deformation of a non-rotating 
star, and treat tab as a flrst order quantity, so that we have: 



STab = SpQab + tab 

The equations we need to solve for equilibrium are thus: 



(2) 



(3) 



where now p ind icates the background de nsity and g{r) = V"!? = GM/r^ . If we expand the perturbation in tensor spherical 
harmonics, as in lUshomirskv et alJ i2000l) . we have: 



tab = trrYlmiraVb - ^Bab) + tr±{r)fab + tA{r)(Aab + ^Yl„^eab) 

where 

eab = Qab - TaTb 

fab = P^^r{ra\7bYlm + rbVaVim) 
Aab = /3" V'VaVbY,,,. + /„6/3"' 

where in equation |H| we have corrected a typo of lUshomirskv et al.l ||200( 
From this expansion one can calculate the quadrupole moment Q22 



Q2 



r 

9{r) 



3 dif'p 
2~dr' 



4 dtr_ 



f3 dr 
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trl 



J Spr'^dr associated with the deformation 



5$ dp 
g(r) dr 



dr 



(4) 

(5) 
(6) 
(7) 
(8) 

(9) 



In lUshomirskv et al.l i2000l) the authors now make the Cowling approximation, thus choose to neglect the perturbations of 
the gravitational potential 5$, and integrate by parts by imposing that the shear modulus p, vanishes at the base of the 
crust. With this condition the surface terms vanish and the final expression for the maximum quadrupole of lUshomirskv et alJ 
J2000l) is 



Q2 



^(4-f/K, + i(6-t/)tA + y|(. 



dr 



(10) 



where U = ^]^^ -I- 2 . It is important to note that we have assumed that the surface terms in the integral vanish because we 
have assumed (as in lUshomirskv et alJ ( 200fi )) that the shear stress vanishes as we approach the base of the crust, however 
this is not an accurate description of the physical situation. The shear modulus in fact increases as one approaches the base 
of the crust from the exterior regions of the star (as can be seen in figure I^J and the transition to = in the core may be 
very sharp; in this case the surface terms may be important. 



2.1 Maximum strain and the Von Mises criterion 

Up to now we have considered the purely elastic response of the crust. However real solids behave elastically only up to a 
maximum strain amax after which they either crack or deform plastically. In order to estimate the maximum size of the 
mountain that can be built we shall assume that the cr ust cracks upon reaching a certain yield strain. Let us investigate this 
in detail, following once again lUshomirskv et al.l (I2OO0I) . We will consider the strain tensor rather than the stress tensor, as 
this allows us to use a simple criterion to establish when the crust breaks. The strain tensor is defined as: 
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0-ab = tab/fJ- (11) 

where fi is the shear modulus. If we define a as — ^(JabO-"*, the Von Mises criterion ijUshomirskv et al.l (|200(t l and references 
therein) then states that the crust will yield when 

> CT-max (12) 

Other criteria exist, such as the Tresca criterion which depends on the difference betwe en the maximum and m inimum 
eigenvalue of the strain tensor, but they are expected to lead to similar qualitative resxi lts j y^h^mir^ky ^t ah* (2000)). If we 
assume the physical variables to be the real part of our complex variables we can write (IUshomirsk^eralT(l200Cl'l ) 

aaba"" = ^a,2.[5R(yi™)]' + a^x[»(/ai,)]' + 4[K(Aa6 + iVimeab)]' (13) 

where arr, <yr± and a a are the components of the strain tensor along the tensor spherical harmonics defined in equation 
By then assuming that the quadrupole is maximised when the equality of equation 11211 is satisfied, i.e. when a = a"max, 
lUshomirskv et alJ (|200(t find that the maximum is obtained when: 

(14) 

0"rX = -y/— ^O'max (15) 

V 5 

CTA = W — g-CTniax. (lOj 

The V on Mises criterion, together with the expression for the maximum quadrupole in (1101 . thus allowed lUshomirskv et all 
i2000l) to conclude that the maximum quadrupole is 

QS'^^ ^ lO^'^g cm2(amax/10-') (17) 

no matter how the strain arises and under the assumption that all the strain is in the Y22 spherical harmonic; strain in other 
harmonics would push the crust closer to the yield point without contributing to the quadrupole. This estimate is thus an 
upper limit for the quadrupole the crust could sustain, and on the energy the star would emit in Gravitational Waves, which 
is related to the quadrupole by: 

. _ 256^Gn« 2 

^gm - -1^-^Q22 (!») 

where Q is the rotation rate of the star. The value (Tmax = 10~^ was chosen becau se it puts Q22 in th e right range to balance 




the accretion torque in the LMXBs, the quadrupole required to do this is in fact JUshomirskv et all llioOC 

Qgn, = 1.6 X 10^*g cm2(300Hz/i.,p,„)''/^ (19) 

if we assume the most basic accretion torque model, as presented for example in I Andersson et "all (I2OO5I) . This would mean 
allowing the spin equilibrium periods of the LMXBs to be set by gravitational waves, and would explain the narrow distribution 
of such spin periods. It is, however, also po ssible to explain the spi n distribution by using a more refined accretion torque 
model, that requires no gravitational waves ijAndersson et all (|2QQ3))- On e should keep in m ind that amax is a very poorly 
constrained parameter, expected to be in the range amax = 10~^ — 10~^ ijRudermanI j 



3 ACCRETED VS. NON-ACCRETED CRUST 

Having reviewed the formalism of lUshomirskv et alJ (l200(j) . let us apply it to the problem of determining the maximum 
quadrupole that an accreted and a non-accreted crust can sustain, and present some of the issues involved. The maximum 
quadrupole defined in equation IIIO^ depends on an integral over the crust involving the quantities g and /i (which is contained 
in the components of the stress tensor tat ~ l-io-ab)- These quantities de pend on t he eq uation of state and composition of the 
crust, which are quite different in the accreted and non-accreted cases llla ensei (l200lh '). The equation o f state, for example, 
determines the thickness of the crust, while the composition of the crust determines the shear modulus iOeata fc Ichimaru.1 
Jl990l)): 



: 0.1194 



(i) (i^~) 

where X„ is the fraction of neutrons outside nuclei, nt is the baryon density, Z is the proton number and A is the atomic 
number. These quantities are plotted in figures and |5| One could then ask which kind of crust could support a larger 
deformation, which is the same as asking which kind of neutron star is potentially a "better" source of gravitational waves. 



4 B. Haskell et al 





Accreted crust 


Non Accreted crust 


Accreted crust 




TOV 


Mass 


IAMq 


IAMq 


1.6Mq 




Radius 


12.56 km 


12.3 km 


12.3 km 




Crust Thickness 


1.76 km 


1.5 km 


1.5 km 




Qmax 


1.6 X 1038 (^) 


g cm^ 2.0 X 1038 (yjpt) g cm^ 


1.1 X 103« (^) 


g cm^ 


e 


1.1x10"'^ 


1.5x10"^ 


7.3x10-* 





NEWTONIAN 



Mass IAMq IAMq 1.7Mq 

Radius 14.1 km 14 km 14.3 km 

Crust Thickness 2.2 km 1.9 km 1.9 km 

Qmax 2.8 X 1038 (^)gcm2 4.2 X lO^* (^) g cm^ 2.5 x lO^* (^) : 

e 1.6x10-'^ 2.4x10"^ 1.1x10"^ 



Table 1. The first two columns show maximum quadrupolc for two stars of equal mass, one with an accreted crust and one with a 
non-accreted crust. The last column is a star with an accreted crust of the same thickness as that of the star with a non-accreted crust, 
but a different mass. Both Ne wtonian equations and the TOV equations are used for the core, in order to illustrate OR effects. The 
equation of state is taken from lDouchin fc Haensell 1200 J) for the core and for the non-accreted crust, while the equation of state for the 
accreted crust is taken from lHaensel &: Zdunikj ll99Cf) . e is the ellipticity ii—iji. 



It has been sugge sted that an accreted crust would sustain a smaller "mountain" as Z is much lower than in a non-accreted 
cr ust JSa tol Jl97 9[)'). To answ er th e question we can evaluate the expression in equation IIUII in the two cases, using the results 
of lPouchin fc Haensell i200ll) and lHaensel fc ZdunikI il990ll for the composition and EOS of the crust. The main obstacle is 
that the composition of the accreted crust is only given up to neutron drip density, because at higher densities the equation 
of state is basically that of a neutron gas, and thus can be m atched to that of a non-a ccreted crust. However, even if the 
EOS is the same, the composition (specifically Z and ^4) is not (iHaensel fc Zduni^ (Il990l) ). In order to obtain an estimate we 
have extrapolated the results available up to neutron drip density to higher densities by noting that the ratio Z/A is roughly 
constant in the two kinds of crust and assuming that it will remain so even in the denser regions at the base of the crust. 
Having done this we can calculate the maximum quadrupole for the two cases. To do this we need to integrate the equations of 
hydrostatic equilibrium for the star (as the core is unperturbed) and we can use Newtonian equations or the TOV equations. 
We shall examine both cases in order to obtain an estimate of how important relativistic effects can be. The question is 
relevant, because even though relativistic effects are not very important in the crust, they do contribute significantly to the 
structure of the star, and thus affect parameters such as the crust thickness which are relevant for our calculation. It would 
then not seem appropriate to use realistic equations of state if we are still describing the core with Newtonian equations of 
equilibrium, as the relativistic effects we are neglecting may be as large as the corrections we expect from the equation of 
state. We will always use Newtonian equations in the crust. 

As we can see from Table [H the non-accreted crust allows for a slightly larger quadrupole, but the results are essentially 
the same for the two models. This is mainly due to the fact that the slight decrease in /x (as seen in figure^ for the accreted 
crust is compensated for by it being thicker. We can try to compare the case of two stars with a crust of the same thickness 
and, as is o bviou s from Table Q in this case there is approximately a factor 2 difference. This is close to what would be 
expected, as lSatol (|l97?i) has predicted that the shear modulus of an accreted crust would be approximately half that of a non 
accreted crust, even if in reality this effect is counterbalanced by the accreted crust being thicker. However it is important to 
note that the result depends on how we extrapolate the composition at the base of the accreted crust. In order to confirm this 
result it would therefore be necessary to have a detailed calculation of the composition of the crust in this region. From Table 
Q it is also clear that relativistic effects can play a role in this calculation, as they appear to make a difference of factors of up 
to a few in the final result (mainly because the crust is thicker in the Newtonian case, and the radius of the star is larger). 



4 VALIDITY OF THE COWLING APPROXIMATION AND SURFACE TERMS 

Following lUshomirskv et al.lll2000ll .in the previous calculation we adopted the Cowling approximation. Such an approximation 
is often used when discussing oscillations, in cases where one expects only the surface layers to take part in the motion. Then 
the high density interior remains essentially unperturbed and one can neglect the effect of perturbations to the gravitational 
potential compared to other forces. We could expect this approximation to be valid also in our (non-oscillatory) context. 
This simplifies the calculation considerably, because it means that the liquid core of the star cannot support any kind of 
perturbation (with no shear and no (5$ there would be no restoring force for a pressure perturbation Sp). This allows us 
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Figure 1. Comparison of the proton number Z and of the ratio Z/A for an a ccreted and a non-accreted crust. EOS from 
iDouchin fc HaenseJ 1200 1^ for the non-accreted case and from lHaensel fc Zdunij ll99(i) for the accreted case. 
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Figure 2. Comparison of the atomic number A and the shear mod ulus /i for an accreted and a non-accreted crust. The EOS is taken 
from Douchin & Hacnscl L2001I) for the non-accreted case and from iHaensel fc Zdunik. 1.1990.) for the accreted case, and then extrapolated 
as described in the text. 



to consider only the crust. In lUshomirskv et alJ l)200dh the authors consider the effect of the approximation and conclude 
that retaining the perturbations of the gravitational potential could lead to changes in the quadrupole f rom 20% to 200%. A 
calculation of the strain that arises in the crust, dropping the Cowling approximation, is carried out bv lCutler et al.l ll2003t) . 
but for the case of a perturbation with a angular dependance, given by rotation. In order to investigate the effects of 
such an approximation for our problem we will present a series of model calculations where the S$ terms have been retained. 
These simple examples can be solved analytically and allow us to estimate how much the Cowling approximation can affect 
the final solution. 



4.1 Constant Density and Shear 

Let us first of all consider a constant density, elastic sphere. In this case we will have 

p{r) = pe{R~r) and n{r) ^ fj,e{R - r) 
H — Const 



(21) 
(22) 
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where Q{R~ r) is the Heaviside step function and we have taken the shear modulus constant over the whole star rather than 
just the crust, in order to make the algebra simpler. R is the stellar radius. In addition to equation ^ we now need to solve 
the perturbed Poisson equation 



(23) 



Let us examine the source term containing 6 p. Take for the components of the strain tensor the maximum from l|16|l . This 
means that the stress tensor is simply tab = fJ-fyab and depends on r only through p,. If we assume that the crust is everywhere 
strained to the maximum and insert our definitions for /i and p into equation ^ we obtain 



dp 



S{R~r)p 
Q{R-r) 

gir) 

S'iR-r) r 



15 

327r. 



167r_ 1 
5 3 

167r 



+ 



167r. 



5{R~r) 



g(r) p\j 5 g{r) 
Returning to Poisson's equation, in the I = m = 2 case we have to solve: 



5$ + 



6 



5$ 



-S{R^ 



rp 



-5^ 



(24) 



(25) 



subject to the conditions that the solution should be regular at the centre of the star, fall off at infinity and be contiimous at 
the surface of the star. The step in the first derivative of 5<E> can be computed by integrating between R — e and R + e and 
then taking the limit e — > 0; this gives 



<5$ (i? + e) - 5$ (R - e) 



(26) 



The solutions of the homogeneous equation are of the form Ar^ + ^ We then clearly have Ar'^ for the interior and B /r^ for 
the exterior. Matching at the surface gives us B = AR^ which can be used to compute the step in the first derivative and 
obtain the interior solution 



5$ ^ ^ Z^*^™" 



/StT 2 



2 i?2p 

We can now substitute this back into equation 12411 and calculate: 



Q22 



47rGp 



R-" 

A-kG 2 p 



sh + I 



(27) 



(28) 



(29) 



which gives us the ratio 
Q.I. ^ 1 

Qsh 2 

where Qah is the part of the quadrupole we would get in the Cowling approximation and Q$ is the part that comes from the 
terms involving 5$, as defined in equation H28|l . It would thus appear that, in this case, by making the Cowling approximation 
we could overestimate Q by a factor of 2. 

Let us see what this means in terms of the ellipticity e = -^7^- We can write the torque acting on a constant density 
ellipsoid as 



32_G 

5c' 



^e\2/5MRyQr 



However we can also write this as 

i 25671- 5 5 5 
J = ——Gil Q22C 
75 

If we compare these two expressions we see that, in the constant density case 
Q22 



= V307r 



3R^M 
4G7rp27?2 



(30) 

(31) 

(32) 
(33) 
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Let us take as canonical values fi = 10'^" dyne/cm^, p = 10^^ g/cmP' , R = IQ km. Substituting these into the above expression 
gives 

-8 



e = 5.3 X 10" 



(fP) (34) 



4.2 Constant Shear in the Crust 

A slight improvement one can make is to impose that the shear modulus vanishes outside the crust. We will still take a 
constant density but we shall take /i of the form: 

At(r) =M[e(i?-r-)-e(r6^r)] (35) 

The calculation is similar to the previous one, but we now have to match the solutions at R and Tb and impose the step in 
the derivative at these two points. It is in fact important to note that as we now have a step at the base of the crust, this will 
lead to surface te rms when we take deriva tives of the shear modulus. These terms are important, as they are not present in 
the calculation of lUshomirskv et alj (l200(jl . where ^ — > at the base of the crust. From this calculation we obtain 

^ lOG^ V~ ^ ' 

^* = "rnc-p V~ ^ ^ 

Again, the ratio between the two is 
Q$ _ 1 

Qsh 2 

and 



(38) 



Q,h = 1.9 X lO^^g cm^ (39) 
In terms of e we have 

which gives us, for the canonical values we have chosen and ri, — 1 km: 

(41) 



/ ^ max \ 



which is nearly an order of magnitude smaller than the previous result for the elastic sphere. So the shell is, as is intuitive, 
"weaker" than the sphere. 

4.3 n=l Polytrope 

A slightly more realistic model is that given by an n = 1 polytrope 

p = Kp' (42) 
which leads to the density profile: 
sin {nr/R) 

p{r) = pc 43 

r-K 

Let us still take ^(r) — /i [0{R — r) — 0{rb — r)] with p constant. The homogeneous Poisson equation thus becomes, in the 
interior: 

5^" + -5$' ~ -^5$ + = (44) 

r r'' R'' 

which admits solutions of the form 

5$ = ^ [{3R^ - ttV^) cos (vrr/i?) + 37r7?r sin (nr/R)] + (45) 

+ ^ [(3i?^-7rV^)sin(7rr-/i?)-37r_Rrcos(7rr/i?)] (46) 

By expanding this solution with a power series around the origin we can see that the regular part is that involving the constant 
B (thus A — 0). This must be matched with the whole solution at rt, and we must then match with the exterior solution 
C/r^ at the stellar surface R. Finally we must impose the steps in the derivative at rt, and R. The resulting solution is quite 
complicated, but leads to 
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^-0.1 (47) 

for a radius _R=10 km and a crust thickness of 1 km. Examining the ratio -S^ for different values of the crust thickness leads 
to very similar results and allows us to conclude that including 5$ always slightly increases the value of the quadrupole. If 
we now calculate an e we obtain: 

^ V30Q ^ V30Q,H ( Q* A , . 

which gives us, with typical parameters for the crust and rj, = 1 km, 

e = 3.1 X 10-' (1^) (49) 
which is an order of magnitude larger than in the constant density cases. 



4.4 Density dependent Shear 

Let us finally take a slightly more realistic model with a density dependent shear. We will take p/ p to be constant which, as 
shown in Figure |3 is not a bad approximation. This gives us 

M = CpQ{r - n) (50) 

where the step function ensures that the shear modulus vanishes in the core {vi, is the radius at the base of the crust). There is 
no need for a step at the surface, as the density goes to zero there anyway. We will take C = 10^® err? j f?" . As in the previous 
section we will take the equation of state to be an n = 1 polytrope, thus obtaining: 

sin(7rr/_R) 

^. = Cpc e(r-ri,) (51) 

The components of the stress tensor will now also depend on r through the shear modulus. 

t„. = ±9(.-n)a^^---^'^(-/^)^^ (52) 
15 r^-K 

4 Cpcsva{-Kr/R)V^R 

tr± = -Q(r~rb)a ' (53) 

5 r^/iY 

4^, X Cpc sin (TTr/R)VSOR 

tx = -er-rtcr-^-^ ^ ' ' (54) 

5 rVTi" 

The procedure we have to follow is identical to that of the previous sections, we have to solve Poisson's equation in the core, 
in the crust and outside the star. We must then impose the boundary conditions at the edge of these regions. These are the 
same as in the preceding sections, i.e. we need to impose the continuity of 5$ and the step in the derivative at the base of 
the crust. 

This allows us to compute the part of the quadrupole that comes from the perturbations of the gravitational potential 
Qif,. For our canonical values we have 

= -2.3 X 10^** g cm^ (55) 

We must also calculate the quadrupole that comes from the shear terms alone. As we have mentioned before this contribution 
involves derivatives of the stress tensor, it also involves delta functions and their derivatives. This means that we have two 
contributions: a surface term that comes from such delta functions and a term that comes from the interior integration. Let 
us examine this case in more detail, as it is a good example of how important the surface terms can be. 
For the stellar parameters given above we find that the integral term is: 

Q>nt = -3.3 X 10^* g cm^ (56) 

while the surface term has opposite sign and gives 

Qe = 4 X 10^* g cm^ (57) 
The two thus add up to 

Qsh = 7.4 X 10^"" g cm^ (58) 

the large cancellations due to the surface term lead to a total quadrupole where the dominant contribution is that given by 
the perturbations of the gravitational potential 
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P ig'cm ) 

Figure 3. The ratio fJ./p, which is taken to be constant in our approximation, is plotted for the case of a non-accreted crust (the case 
of an accreted crust is qualitatively similar). As we can see it is approximately constant to within factors of two. 



Qt. 



-1.5 X lO''** g cm^ 



(60) 



corresponding to an e of 



Which is smaUer than the previous resuh and is a direct consequence of the canceUations that occur at the base of the crust. It 
is thus quite clear that the correct treatment of the boundary can have an important eflect on this calculation. In particular it 
would appear that one cannot take /i — > at the base of the crust, as the step is actually quite sh arp (as shown in Figure 
and we have seen that this can lead to large surface terms. We can, however, conclude, similarly to lUshomirskv et al.l (jEoOOl ). 
that including the perturbations of the gravitational potential 5^ will affect the results by a factor 0.5-3. Thus the Cowling 
approximation is not consistent, as the corrections given by Sij) are not negligible. It must, however, be noted that when we 
drop the Cowling approximation we are no longer guaranteed that the components of the strain tensor calculated in Q will 
give us the maximum quadrupole. 



5 PERTURBATION FORMALISM 

Up to now we have followed the method of lUshomirskv et alJ i2000l) and simply taken the crust to be strained to the maximum 
throughout. However, as we have seen previously, the boundaries can have a strong impact on the result, so in order to carefully 
impose the boundary conditions at the surface of the star and at the base of the crust let us reformulate the problem in terms 
of perturbations in the crust. We shall impose (in a si milar way as w as done to calculate rotational deformations in the 
slow-rotation approximation, as presented for example in iTassouH l)l978l) 'l that the star has a small quadrupolar deformation, 
such that we can define a new variable a which is a function of the radial variable of the spherical background model 

a{r)^r[l + ^{r)Y22ie,(l>)]. (62) 

We shall assume that a(r) labels the equipotential surfaces of the gravitational potential (which will thus have the same 
deformed shape). It is important to note that in the core these surfaces also describe the "isobaric" and "isopycnic" surfaces. 
This is a consequence of the hydrostatic equations. In the crust, where we also have elastic terms, this is no longer true. In the 
crust, surfaces of constant $ are NOT surfaces of constant density, so while e is a good variable to use for the integration, it 
does not represent the "shape" of the star, which can be obtained by calculating the perturbations of the density (Sp) or of the 
pressure (as we will assume the equation of state to be barotropic). We expand the displacement vector along the orthonormal 
basis formed by the radial unit vector f, ^VYim and x VYim- This allows us to use a fiat background 3-metric gat- We 
shall ignore the axial part, proportional to f x VYim as this does not couple to the scalar perturbations such as Sp which are 
proportional to Yim- The displacement vector thus takes the form: 

r = er(r)Y,,„r" + a VY,™ (63) 
If we assume that the crust responds elastically it is useful to define the stress-energy tensor of the solid Tik by 

Tik = —pgab + tab (64) 
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where we have defined the trace-free tensor 

tab = M {^aib + '^bia ^ ^ffab V^c^ (65) 

which describes the shear stress in the solid. This allows us to write the equations of static balance as 

V'tab = pV6$ + VbP (66) 

which we can also write in the form 

V"t,6 = pV5($ + H) (67) 
where we have introduced the enthalpy 

"-Im 

We also have to supply a barotropic equation of state 

P = P{p) (69) 
The remaining equations we need are the continuity equation and Poisson's equation 

&P = -drpir- pdrir+p-U-'i-ir (70) 

r r 

= 4ttGp (71) 

and the condition that the unit vector a'' given by the new radial variable is orthogonal to the equipotential surfaces of $, or 
parallel the gradient of "I>, i.e. 

Va X V* = (72) 

It is important to note again that a is not related to the 'shape' of the star (as the deformations of the density are) because 
we now have elastic terms in our equations of force balance. However the introduction of the variable ^ makes the equations 
for 5<E> more tractable numerically. We will now carry out first order perturbations, treating as a first order quantity. This 
corresponds to assuming that the background spherical configuration is at zero strain. We could imagine this to be plausible 
for a star which has been spun down by dipole radiation and been in an essentially non rotating configuration long enough 
for plastic fiow to relax the crust. If this star then starts accreting it would be spun up and correspond to the situation we 
are describing. It is of course not relevant for the case of a spinning down star with a crust that formed while it was rapidly 
rotating and for which the reference zero-strain configuration would not be spherical. This problem is, however, a much more 
complicated one, as it would require us to know the reference shape of the star, which will in general depend on the history 
of the crust. 

We shall assume that every perturbation has a Y22 angular dependance, i.e. is of the form 5p{r, 0, (/>) — 5p{r)Y22{Q, <!>)■ This 
is because we are interested in the quadrupolar part of Sp which will give the dominant contribution for gravitational waves. 
However it is worthwhile noting that, even though they do not contribute to the gravitational wave emission, perturbations 
with a different angular dependance (e.g. a 120 deformation due to rotation) still build up stress in the crust, and bring 
it closer to the breaking point, thus allowing a smaller quadrupolar deformation. For example, there could be a strong Y20 
deformation, as the star is rotating, and this could build up high levels of strain in the crust. It is not however easy to evaluate 
this quantitatively, as the exact level of strain depends on the reference shape of the crust (the unstrained configuration), 
which in turn depends on the history of the star. Let us then write the perturbation equations we need (/3 = ^yT(T^-T) = y/6): 

Vbt"^ = pV" [{5H + (5$)F22] Hydrostatic equifibrium 

5p — —drP^r — pdr^r + p^^± — 2^^,- Continuity 

5p = ^^^j^Sp = c^Sp Equation of state 

SH ^ ^ Enthalpy 

(5<1> — —tprdr^ Level surfaces (from eg. 17211 ) 



Poisson 
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This gives us 7 equations for our 7 variables ^r, ^x, Sp, Sp, SH, 5$. Some of these equations are simply algebraic relations 
defining some quantities (such as the enthalpy). We can thus re- write the equations as three second order equations containing 
the variables ^r, ^± and ip. Once we have solved for these variables we can compute Sp, for example from the continuity equation. 
Let us write the stress tensor tat in terms of our displacements: 



dr r 



(-2/3^ + 2/3^ + 2/3^) + fah (: 



dr 



dr P P 

Working out some of the algebra we obtain for Vbt'^'', breaking it down in its two components along f" and V''l22 
dfi f d^r 2 



Ta^bt 

and 

Val22V6t'" 



Ylrn 

3r2 



'fr i^'^^ + '^/^W2-re.)+M 



dP£,r 2 



^Pr + 7Cx/3 + 8^r - 8^,- - 
dr dr 



1 

/3r 



dr 



{rPir + r'^+ru)+^,{^ 



^^r/3/3-4/3/3^ex+8/3C./? + 2^. + .^^, 



The hydrostatic equilibrium equations then give 



fa^bt" 



pdr{SH + S^) 



(73) 
(74) 

(75) 
(76) 



(77) 
(78) 



Note that the term involving 5pdr{H + $) is not present because the background quantity dr{H + $) vanishes in equilibrium 
(as the background spherical configuration is a zero-strain one). At this point we can write the equations as a set of four 
algebraic equations and three coupled second order equations for the displacements: 



^ (3c?p + 4m) pr' = - (2p + 7p^^P + iclpp" - 3^c?) 



(79) 



.dp o «2 , g 2 d^p 2 o 2 (dpy 2 , r,d 2 2 r. 2 2 \ c 

4:p—r-8pp-3ppP +3c^p-^r " ( — I r +b — c^rp -6c^p j 



d^p 



dp^ 



d^r 

dr 



- (sppr + Sc^p^r^ + ip-j^r^ + 3 — c^r^p 
V dr dr dr 

+ {ppPr + SciPp\) ^ + (3r^^^p) ^ - (sr^A + Ur'p^c) ^ 
^ ' dr \ dr dr J dr \ dr J 



d^ 

dr^ 



[pr^] = - (-4^/3' + 3^ - 3c^/3'p) U - {^'^^P + 8^/? + 3/3c'^r + 6c^/3p) 



dp 



I dp 



dp 



dii 



>.dp 



'^r^ i'^'s'i) = i^-GPp^) ex - - 8.Gp^) + {^r^cf^P ^ P^4 + S.GpV) ^ 



'.dp'' 



dp 



dp 



(80) 



(81) 



dr 



If we make the Cowling approximation these equations reduce to those used bv lMcDermott et al.l 



5.1 Boundary conditions 

Let us now move on to the key issue that was ignored in our previous analysis. We need to prescribe boundary conditions 
at the surface and at the base of the crust. To do this in a consistent manner we should solve the perturbation problem in 
the core and then impose the continuity of the gravitational potential and the continuity of the perpendicular and radial 
components of the traction vector. We assume that the perpendicular tractions vanish outside the star and also in the liquid 
core, as it cannot support shearing. The core could in principle sustain pressure and density perturbations, and these will 
appear in the radial component of the traction. The definition of the traction vector is 

t" = r^'db (82) 
where T"* is now the complete stress tensor, i.e. 

T"' = -Spg"'' + t"'' (83) 
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This gives us 



e = r 



(^dir \ 2 (dp dir , yja^ , ( . ^ , 



(84) 



Imposing that the component of this along ^^22 must vanish gives us the condition (at the base of the crust): 



^(,,) = _M!:^_^ (85) 

dr rb rb 

We then assume that the radial component of the traction is continous, thus equal to tcore- This leads to 

tcore = ['^^(rb) - + -^f3 j + c.(r,) (^^(r,)C.(r,) + p— (r,) - j (86) 

We shall consider the full perturbation problem in the following sections, but first, for simplicity, let us consider some 
simple boundary conditions: we shall assume that the fluid core is unperturbed and remains spherical. Such a condition is 
not physical, but will allow us to consider a full GR fluid core (as we only need the background model) and quantify how 
important relativistic effects can be. However, it is not a natural condition since it needs a force acting at the surface as well 
as some forces to be acting at the base of the crust to keep it spherical. This is why we cannot impose the traction conditions 
but will rather impose that the radial component of the displacement (fr) vanishes (i.e. the core remains spherical and does 
not 'move' into the crust) and that the shape of the inner edge of the crust remains spherical, which implies that dp must 
vanish (remember that the background is already spherical, so there must be no 122 perturbations if it is to remain spherical). 
We also impose continuity of the gravitational potential, and thus impose e{rb) = (which corresponds to 5$ vanishing at 
the base of the crust). Our boundary conditions are thus 

& = (87) 

= ^ii_^,^^_2ii (88) 
dr rb dr rb 

■0 = (89) 

The shape can in fact be described by a variable deflned to label the equipotential surface of the density or pressure (very 
much like tp)- We can thus define 

^ - (90) 



rdrP 

At the surface we may impose that the shape is given by a certain Sa = es(-R). Wc will then search for the value of Es that 
causes the crust to crack (at some interior point, not only at the surface), following the Von Mises criterion. The final thing we 
need, before considering the surface of the star, are the background quantities for our model, the density and the derivatives 
of the gravitational potential. 

These calculations are all Newtonian, and thus not very accurate in the core where we expect relativistic effects to be 
important. As we are imposing that the core is not perturbed we can integrate the TOV equations in the core, and then 
match them to the Newtonian quantities wc use in the crust (at the low densities of the crust the Newtonian approximation 
is reasonably accurate). This is important as it will allow us to use a realistic equation of state for the core, which would be 
pointless if we were considering a fully Newtonian star. It is in fact well known that stellar models with the same EOS and 
central density can give very different results in the Newtonian and relativistic cases. 

As is obvious from the above equations, both the background quantities p and p vanish at the surface of the star, thus 
causing the tractions to trivially vanish. We shall then need to impose regularity conditions. Starting from our background 
model we can expand our equations around r = R (the surface) and thus obtain 



dr'^ (r - R) 

d^a 1 



_ _der 

R 2 R dr 
R^'^ R^ dr 



+ 0{l)+0{r-R) + .. (91) 
+ 0{l) + 0(r-R) + .. (92) 



dr2 (r - R) 

We now have two second order equations for and ^± which have decoupled from the equation for ijj (which wc are not 
interested in, as we already know what boundary condition to use for ^). The solution involves a combination of Bessel 
functions, this allows us to discard the ones that diverge at r = J? (the Bessel functions of the second kind) and thus be left 
with a solution described by two constants, A\ and A2. Expanding around r = we have: 

= fi{Ai,A2)+gi{A^,A2){r-R) + 0{{r-Rf) (93) 
a = h(Ai,A2)+g2{AuA2)(r-R) + 0{(r-Rf) (94) 
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This allows us to impose relations between the functions and their derivatives close to the surface: 

which is the boundary condition we impose close to the surface, together with the magnitude of the perturbation, i.e. the 
value £3. We thus have a total of three conditions at the surface and three at the base of the crust. 

5.2 Perturbations of the core 

So far we have derived boundary conditions for the problem of an unperturbed core with a perturbed crust. We shall now 
turn our attention to the more physical problem of a perturbed crust with a perturbed core. Let us solve the whole set of 
equations in the core and the crust. In the core we need to solve the same set of equations as in the crust, but now we have 
no shear modulus, and so they reduce to simple hydrostatic equations 

Vp = -pV$ (97) 

which we can expand along our basis vectors 

r'y;™ {drSp + pdrS^ + 5pdr^) = (98) 
VY;™ {5p + p<5$) = (99) 

The equations simplify to 

5p = -p5$ (100) 
Sp = (101) 
We can substitute these into Poisson's equation 

V^^* = AttGSp (102) 
and obtain an ODE for 5$ 

—— + -- ^5$^47rG-^5$ 103 

which we can solve together with equations (IIOH . The only boundary conditions we need to impose are those of regularity at 
the centre, i.e. that the solution close to r = is of the form Ar^ where A is a constant. This leads to 

d(5$ (5<1> 

= 2 — close to the centre (104) 

dr r 

We then impose continuity of the tractions at the base of the crust. 



6 RESULTS FOR AN N=l POLYTROPE 

Let us consider a specific stellar model for which we take an n = 1 polytropic equation of state p — Kp^, which leads to 

5p = 2Kp5p (105) 

As in previous examples the shear modulus is taken to be such that fi/p — 10^^ err? j . We solve the equations numerically by 
using a shooting method with a Runge-Kutta 5th order variable step size routine (Cash-Karp). We calculate the quadrupole 
moment and gradually increase the value of t/) until the crust breaks (following the Von Mises criterion). The equations are 
integrated in three diff'erent cases: 

(i) Newtonian equations throughout the star, only perturbations of the crust 

(ii) TOV equations in the core, Newtonian in the crust, only perturbations of the crust 

(iii) Newtonian equations throughout the star, perturbations of the core and crust 

In the first two cases the core is unperturbed, so we impose the boundary conditions derived in section 5.1, the last case includes 
a perturbed core, so the boundary conditions are the ones described in section 5.2. We choose these cases as an unperturbed 
core has the advantage that it can be obtained by integrating either the Newtonian equations of hydrostatic equilibrium, or 
the TOV equations. This allows us to estimate how important overall relativistic effects may b e. The results are presen ted in 
Figure|l|taking different densities for the base of the crust. This allows us both to compare with lUshomirskv et alJ i2000l) . who 
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6.8xl03''(e 


= 1.4 X 10" 
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Table 2. Maximum quadrupole, in tlic case of (Tmax = 

10-2 for an n=l polytrope, and the base of the crust at p^ = 1.6 X 10^'^ g/cm^ 
and p h = 2.1 X lO^** g/cm^. The parameters of the stellar models are given in table IHlfUCB indicates the result of lUshomirskv et alj 
J2000|) V 



Mcore (Mq) crust thickness (km) K (g -'-cm^s ^) pc (g/cm^) p;, (g/cm"^) 



Tov^ 


-Newt 


1.358 


1.115 


7.24x10-* 


2.45x10*-'^ 


1.6x10^5 


Newt 




1.37 


0.68 


4.25x10** 


2.2x10*^ 


1.6x10*5 


TOV^ 


-Newt 


1.33 


1.39 


7.07x10* 


2.54x10*^ 


2.1x10*5 


Newt 




1.35 


0.88 


4.25x10* 


2.2x10*^ 


2.1x10*5 



Table 3. Parameters of the two stellar models, the fully Newtonian model and the relativistic model with a Newtonian crust, for a 
density at the base of the crust of p = 1.6 X 10** g/cm^ and p = 2.1 X 10** g/cm^ 



take t he density at the base of the crust to be 2.1 x 10**g/cm^, and to use the value 1.6 x 10*'*g/cm'', suggested bv iHaense] 
i200ll) . The results for these two particular cases are shown in Table |5| The Newtonian calculation with the unperturbed core 
can be fitted to the formula 

The dependence on pt of the calculation with the relativistic core is not well fitted by a simple power law, as we can see from 
figure |1] where it is evident that for thicker crusts the maximum quadrupole is growing at a slower rate. However, if we take 
Pb ~ 1.6 x 10** and R — 10 km, the dependence on the stellar mass is well fitted by 

/ \ -0.3 

Omax = 1.6 x lO'" { j g cm^ (107) 

The Newtonian calculation including the perturbations of the core is also very interesting. As we can see from Figure |1] the 
quadrupole built up in the core is predominant if we take the base of the crust at low densities, which is equivalent to having 
a thin crust. The crust quadrupole, on the other hand, becomes predominant for thick crusts, and is well approximated by 
the Newtonian calculation with the unperturbed crust, as we can see in Figures 0] and |^ It is also very instructive to observe, 
in the case of the perturbed core, how the maximum quadrupole depends on the mass of the star (Figure (HJ. Unlike the other 
two cases the quadrupole gets larger as the star gets heavier. This is because the quadrupole from the perturbations of the 
core is growing as the star is made more massive, while the crustal quadrupole is in fact decreasing. The dependence on the 
stellar radius is shown in Figure El for the case of a perturbed Newtonian core and that of an unperturbed relativistic core; 
the two cases are similar and not altogether surprising. All calculations are linear in the breaking strain. 

It would appear from the results of our calculations that considering the perturbations of the core can make a difference 
in the result, mainly for thin crusts. On the other hand it appears that if the crust is thick, then the crustal quadrupole 
plays a predominant role, and is well approximated by taking a perturbed crust on top of an unperturbed core. In this case 
relativistic effects can become important, and it would then not seem appropriate to use realistic equations of state, but still 
describe the core with Newtonian equations of equilibrium. On the other hand it would seem to be the case that one can 
consider an unperturbed relativistic core (with a realistic EOS), then use a realistic EOS and p in the Newtonian crust and 
still have a reasonable approximation to the full problem. We are of course still solving the problem in Newtonian elasticity, 
so not solving the problem in "full" GR, but at the densities of the crust the Newtonian approximation is a reasonable one, 
and makes the problem considerably simpler than it would be in relativistic elasticity. 



6.1 Application to the Accreted vs. Non-Accreted case 

We can now consider the stellar models presented in Tabled for two stars, one with a non-accreted crust and one with an 
accreted crust. We shall first take the two stars to have the same mass,then we shall also consider two stars with the same 
crust thickness. As the crust is now thick, as shown in Tabled we expect the crustal quadrupole to dominate and can thus 
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Figure 4. On the left hand side we have the maximum quadrupo le plotted first with respe ct to the density of the base of the crust. We 
present the three cases described in the text and also the result of lUshomirskv et al.l 1200(1) . marked as UCB. On the right hand side we 
have the maximum quadrupole, for the three cases considered in the main text, plotted with respect to the crust thickness 




consider the approx i matio n of an unperturbed core to be reasonable. This then allows us to use the equation of state given by 
IPouchin fc Haensell ll200lll and solve the TOV equations in the core and the Newtonian equations of hydrostatic equilibrium 
in the crust, thus taking relativistic effects into account. The results are shown in Table 01 As expected the values for the 
quadrupole do not differ significantly. We also confirm that a neutron star w ith an non a ccreted crust could provide a slightly 
better source of gravitational waves, as was expected from the literature iSatd il97Sl) '). The values of the quadrupole (or 
equivalently the ellipticity) are now larger than in the previous case as we are taking a realistic model for the shear modulus 
(as shown in figurej^J. This leads to values of ^ which are significantly larger than in the constant ratio /j,/ p approximation, 
especially in the higher density regions close to the base of the crust, which give rise to most of the quadrupole. 

The values we find for e are now comparable with the upper limits of e « 10~^ set by LIGO for some of the closer pulsars 
iAbbot et ai \200^ ). The best limits, however, are still set by pulsar spindown, and can be seen in Figure|7| compared to the 
limits set by this work. 
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Figure 6. Dependence of the maximum quadrupole on tlie stellar radius, in the case of a perturbed Newtonian core and of an unperturbed 
TOV core. 



Accreted crust 



Non Accreted crust 



Mass 
Radius 

Crust Thickness 
Qmax 



1.4M0 

12.56km 

1.76km 

1.8 X 10^9 (- 

1.3x10"^ 



— jg cm^ 



1.4M0 
12.3km 
1.5km 
3.1 X 10^ 
2.4x10"' 



Accreted crust 



1.6M0 
12.3km 
1.5km 

1.6 x 10^8 (- 
l.lxlO"'^ 



Table 4. Maximum quadrupole for two stars of equal mass, one with an accreted crust and one with a non-accreted crust. The last 
column shows a star with an accreted crust of the same thickness as that of the non accreted c rust, but a different mass. T he perturbative 
formalism of section 5 is used and the base of the crust is taken to be at 1.3 10^*g/cm'^ as in lDouchin fc: HaenseJ 120 



7 CONCLUSIONS 

We have reviewed the method of lUshomirskv et al.l (I2OO0I1 to determine the maximum deformation that a Neutron star crust 
can sustain. This procedure has then been apphed to the case of an accreted and a non-accreted crust, in order to estimate 
whether one of the two would produce a better source of gravitational waves. It appears that the two kinds of crust are 
essentially equival ent, with th e non-accreted crust being able to support a slightly larger quadrupole, as has been suggested in 
the literature (^e.g. lSatol l)l979l) '). It is important to note, however, that to obtain this result is has been necessary to extrapolate 
the composition of the crust in the high density regions at the base of the accreted crust. In order to confirm this result it 
would therefore be necessary to have a precise model for the composition in this region. 

We have then examined a series of toy models, to determine what effect the Cowling approximation can have on the 
problem and, more crucially, if it is safe to assume that the shear modulus /i — > at the base of the crust. We conclude that 
this is not the case, as the shear modulus is expected to vanish over a thin transition region close to the base of the crust, 
which can lead to large boundary terms. The Cowling approximation can also make a difference of factors of a few in the final 
result. 

In order to drop the Cowling approximation, and have better control over the boundary conditions, we have developed 
a perturbation formalism and applied it to the case of an n = 1 polytrope with a co nstant ratio ^/ p. The results of these 
tests have confirmed that the Cowling approximation used in lUshomirskv et al.l J2000l) can have a large impact. Our results, 
however, also suggest that it is not a bad approximation to consider an unperturbed core and perturb the crust, if the crust 
is sufficiently thick. This allows us to solve relativistic equations in the core, which can give corrections that are significant if 
we intend to use a realistic equation of state for the star. Having developed the formalism we apply it to the original problem 
of the accreted and non-accreted crusts, with a realistic model for p, thus confirming that the two crusts can sustain very 
similar quadrupoles, the non accreted crust being able to sustain a slightly larger "mountain". What is surprising is the fact 
that the quadru pole found with this form alism is approximately an order of magnitude greater than the results obtained with 
the formalism of lUshomirskv et al.l (l200fi) . Such a significant quadrupole, if it were confirmed, may allow us to put an upper 
limit on the breaking strain of the crust from observations. The problem, however, is not quite so simple. The calculations we 
have presented give a "maximum" quadrupole in as much as we have assumed all the strain to be in the Y22 harmonic. But 
the stars we consider are rotating and must therefore be rotationally deformed, thus putting strain also in the 1^20 harmonic. 
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Figure 7. The upper limits on the eUipticity e set by pulsar spin down (taken from the ATNF pulsar database, 
www.atnf.csiro.au/ research/pulsar/psrcat/) and the thoretical upper limits set by the perturbative method this work (top line) and 
with the method of lUshomirskv et alj 120001) (bottom line). 



as spin-up or spin-down moves the star away from its reference shape. This could bring the crust much closer to the yield 
point and significantly decrease the size of any additional "mountain" that it could sustain. Unfortunately, estimating exactly 
how much strain is built up by rotation is a complicated problem, as the unstrained configuration of the crust depends on 
the history of the crust, and there may have been many starquakes during the spin-down and subsequent spin-up phases. 
However, if a neutron star can sustain a quadrupole this large for a breaking strain of crmax = 10~^, it would appear that 
even if the breaking strain is smaller there is still a possibility that the quadrupole may be large enough for gravitational 
waves to have a role in the spin evolution of the star. 

Finally we should note that our approach was chosen in order to obtain a "maximization" argument, independent of the 
mechanism that produces the deformation, and to show that o verall relativistic effects can be important. We have, therefore, 
integrated a set of source free equations, analogous to those of lMcDermott et al.l lIlQSSi) . with the deformation as a boundary 
condition. It would, however, be important to understand the physical processes that caus e the deformation and int egrate the 
same equations with a source, such as the temperature or composition variations used in lUshomirskv et alJ ll2000l) . Another 
important effect we are neglecting is that of the magnetic field. It is well known that the equilibrium shape of a star will not be 
spherical in the presence of a magnetic field and, even though crustal effects may be dominant early on during accretion, the 
process of magnetic burial cou ld make the field locally s trong enough for the magnetic deformations to become important. This 
effect was recently studied bv lMelatos fc Pavnj liooi) and lPavne fc Melatoj 1120051) ■ who were the first to perform globally 
self-consistent MHD calculations of magnetic mountains on accreting NS as gravitational wave sources. These calculations 
did not include crustal effects (as this would make the problem much less tractable), nevertheless it is important to note here 
that if one were to solve the full inhomogeneous problem, including the magnetic source terms (or terms due to other physical 
process we have neglected), one could obtain deformations larger than the crustal maximum we have calculated. 

It would also be of some int erest to consider the full problem in a relativistic formulation of elasticity, such as that 
of iKarlovini fc SamuelssonI i2003l) . as we have shown that the perturbations of the core can be important, especially in 
determining the dependence of the maximum quadrupole on the stellar parameters. We intend to turn our attention to these 
problems in the near future. 
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